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ABSTRACT 

The observed relationship between X-ray luminosity and temperature of the diffuse intercluster medium 
clearly shows the effect of nongravitational heating on the formation of galaxy clusters. Quasar feedback into 
the intergalactic medium can potentially be an important source of heating, and can have significant impact on 
structure formation. This feedback process is a source of thermal Sunyaev-Zel'dovich distortions of the cosmic 
microwave background. Using a simple one-dimensional Sedov-Taylor model of energy outflow, we calculate 
the angular power spectrum of the temperature distortion, which has an amplitude on the order of one micro- 
Kelvin. This signal will be at the noise limit of upcoming arcminute-scale microwave background experiments, 
including the Atacama Cosmology Telescope and the South Pole Telescope, but will be directly detectable with 
deep exposures by the Atacama Large Millimeter Array or by stacking many microwave images. 
Subject headings: cosmology — S-Z effect: AGNs — CMB 



1. INTRODUCTION 

The tiny temperature fluctuations in the cosmic microwave 
background, arising from small density perturbations in the 
early universe, have proven to be our most powerful con- 
straint on properties of the universe (Spergel et al. 2006), 
spectacularly fulfilling theoretical expectations (Jungman et 
al. 1996a). The WMAP satellite, the culmination of nu- 
merous microwave experiments, has mapped the microwave 
sky with an angular resolution of 15 arcminutes (Hinshaw et 
al. 2006, Jarosik et al. 2006). Attention is now turning to 
smaller angular scales of a few arcminutes, where temper- 
ature fluctuations arise due to interaction of the microwave 
photons with matter in the low-redshift universe (for a review, 
see Kosowsky 2003). These secondary anisotropics arise 
due to nonlinear evolution of gravitational potentials (Rees 
& Sciama 1968), peculiar velocities of ionized dense regions 
(Ostriker & Vishniac 1986), gravitational lensing (Blanchard 
& Schneider 1987, Bernardeau 1997), and microwave scatter- 
ing from hot electrons (Sunyaev & Zeldovich 1972). The last 
of these, now known as the Sunyaev-Zeldovich (SZ) Effect, 
provides a powerful method for detecting accumulations of 
hot gas. 

Galaxy clusters, which contain the majority of the thermal 
energy in the universe, provide the largest SZ signal, and up- 
coming microwave surveys like ACT, the Atacama Cosmol- 
ogy Telescope (Kosowsky et al. 2006), and SPT, the South 
Pole Telescope (Ruhl et al. 2004), are expected to detect 
thousands of clusters over their entire range of redshift, po- 
tentially placing strong constraints on the evolution of struc- 
ture in the universe (e.g., White, Efstathiou, & Frenk 1993, 
Bahcall, Fan, & Cen 1997, Molnar et al. 2004, but see also 
Francis, Bean, and Kosowsky 2005). A number of other as- 
trophysical processes will also create measurable SZ distor- 
tions, including peculiar velocities during reionization (Mc- 
Quinn et al. 2005, Iliev et al. 2006), supernova-driven galactic 
winds (Majumdar, Nath, & Chiba 2001), black hole seeded 
proto-galaxies (Aghanim, Ballad & Silk 2000), kinetic SZ 
from Lyman Break Galaxy outflow (Babich & Loeb 2007), 
effervescent heating in groups and clusters of galaxies (Roy- 
chowdhury, Ruszkowski & Nath 2005) and supernovae from 
the first generation of stars (Oh, Cooray, & Kamionkowski 
2003). Here we investigate one other generic source for SZ 



signals: the hot bubble surrounding a quasar or active galactic 
nucleus. Previous related studies were done by several au- 
thors (e.g., Natarajan & Sigurdsson 1999, Yamada, Sugiyama 
& Silk 1999; Lapi, Cavaliere, & De Zotti 2003; Platania et al. 
2002). 

Analytic studies and numerical simulations of cluster for- 
mation indicate that the temperature and the X-ray luminos- 
ity relation should be L x ~ T 2 in the absence of gas cooling 
and heating (Peterson & Fabian 2006). Observations show in- 
stead that L x ~ T 3 over the temperature range 2-8 kev with a 
wide dispersion at lower temperature, and a possible flatten- 
ing above (Markevitch 1998, Arnaud & Evrard 1999, Peter- 
son & Fabian 2006). The simplest explanation for this result 
is that the gas had an additional heating of 2-3 keV per par- 
ticle (Wu, Fabian & Nulsen 2000; Voit et al. 2003). Several 
nonthermal heating sources have been discussed in this con- 
text (Peterson & Fabian 2006), and quasar outflow (Nath & 
Roychowdhury 2002) is one possibility. The effect of this 
feedback mechanism on different scales of structure forma- 
tion have been addressed by several authors (e.g., Mo & Mao 
2002; Oh & Benson 2003; Granato et al. 2003). Recent 
progress in understanding the relation of black hole masses 
to dark matter haloes (Merritt & Ferrarese 2001; Tremaine 
et al. 2002) provides a simple connection between individual 
quasars and their overall cosmological impact; Scannapieco 
& Oh (2004) have studied the global impact of quasar feed- 
back into the IGM on structure formation in the standard cos- 
mological model. In this paper, we use this model to inves- 
tigate the quasar feedback process as a source of microwave 
background SZ distortion, and discuss its detectability by up- 
coming SZ surveys and millimeter wave experiments. 

2. QUASAR OUTFLOW MODEL 

A quasar injects a substantial amount of energy into the sur- 
rounding gas while it is active. Following Scannapieco & 
Oh (2004), we assume the black hole powering the quasar 
shines at its Eddington luminosity and returns around 5% of 
this energy to the galactic gas, eventually disrupting its own 
fuel source after a dynamical time of the cold gas surround- 
ing the black hole, fd yn ~ 5 x 10 7 (1 +z)" 3//2 yr (Barkana & 
Loeb 2001). The efficiency factor is estimated by requiring 
that the energy is sufficient to explain the preheating inferred 
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from galaxy cluster X-ray observations. The duration of the 
energy injection fd yn is much shorter than the expansion time 
of the resulting bubble of hot gas (on the order of 10 9 years), 
so we assume an instantaneous point source injection of en- 
ergy into the intergalactic medium. The total energy input is 
just the product of the luminosity, the efficiency factor, and 
the duration. The dependence of the luminosity on black hole 
mass can be converted to a dependency on galaxy halo mass 
via the M^-a relation (Merritt & Ferrarese 2001; Tremaine 
et al. 2002) and a relation between galaxy rotational velocity 
and velocity dispersion (Ferrarese 2002). The resulting total 
energy injected for a quasar turning on at redshift z in a halo 
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of mass Mhaio is approximately E = O.ObM^ (1 + z), where 
Mi 2 =M halo /lO 12 M . 

For simplicity, we assume that after the energy injection, a 
hot bubble evolves adiabatically and expands into a medium 
of uniform overdensity. The one-dimensional Sedov-Taylor 
solution can be used to model the radius and temperature of 
the region contained by the blast wave (Scannapieco & Oh 
2004). The radius scales as 
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where E^o is the energy in the hot medium E in units of 10 60 
ergs, 5 S is the ratio of the density of the surrounding medium 
to the mean cosmic baryon density, and foyr is the expansion 
time of the bubble in units of 10 9 years. We have assumed 5 S 
equal to one which means, an ambient baryon density equal to 
the cosmological density. Many quasars may reside in denser 
environments, but the y-distortion scales only very weakly 
with density. We also have assumed an exact proportionality 
between halo and black hole mass; some quasars may have 
substantially higher energy outputs. The temperature of the 
bubble scales as 
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The density of the gas inside the bubble is assumed to be 
uniform and equal to the density of the gas outside the bubble. 
The actual density profile varies with radius (e.g., Shu 1992), 
but not strongly, and for simplicity we assumed it to be con- 
stant. The cooling time of the gas is of the order of a Hubble 
time, so we ignore radiative cooling (cf. Figure 2 of Scanna- 
pieco & Oh 2004). Energy loss into compression against gas 
pressure and gravitational potential energy changes are also 
ignored. We also assume that all quasars eject their energy at 
a single redshift, z;„. A more realistic model would integrate 
over a range of input redshifts, but this simple assumption will 
give the correct order of magnitude for the resulting SZ distor- 
tions with this mean redshift, since the SZ effect is essentially 
independent of redshift. 

3. THE SUNYAEV-ZELDOVICH DISTORTION 

The Compton y-parameter characterizing the SZ spectral 
distortion is proportional to the line-of-sight integral of the 
electron pressure. Since each individual source is assumed to 
be spherically symmetric, with constant temperature T e and 
number density n e of electrons inside the hot bubble surround- 
ing the quasar, the y-distortion y(9) on the sky will be az- 
imuthally symmetric, depending only on the angle 9 between 
the bubble center and a particular line of sight. The angular 
Fourier transform of the y-distortion is given by (cf. Peebles 
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with oj the Thompson cross section and £>a(z) the angular 
diameter distance to redshift z. The square-root term in the 
integrand is the path length through the hot bubble at an angle 
9 from its center. This integral can be performed analytically 
(Gradshteyn & Ryzhik 1980) to give 
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Note that T e and R s depend on both the halo mass M and the 
redshift z, and n e depends on z. 

The y-distortion on the sky can be conventionally expanded 
in terms of the spherical harmonics as y(n) = Ydm a lmYlm^)- 
The angular power spectrum is then obtained as C\ — ( | a/ m | ) , 
an ensemble average over the coefficients. The power spec- 
trum has two components (Cole & Kaiser 1988, Cooray & 
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where Cf is the contribution 



from Poisson noise of the random galaxy distribution, and Cj 
comes from the correlation between galaxies. The two terms 
are given as (Komatsu & Kityama 1999; Majumder, Nath & 
Chiba 2000) 
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where 

HM,z) = dn{ ^' n) b{M,z m ) yi {M,z), (7) 

k[(z) = 1/Da(z) is the wavenumber corresponding to the mul- 
tipole angular scale I at redshift z, dV /dz is the comoving 
volume element, dn(M,z)/dM is the differential mass func- 
tion, P m (k,z) is the matter power spectrum, and b(M,z) is the 
linear bias factor. The expression for the correlated piece uses 
the Limber approximation. 

To determine the cosmological functions, we assume a stan- 
dard ACDM cosmology with Q m = 0.31, Q h = 0.044, and a 
Harrison-Zeldovich primordial power spectrum P(k) = k. The 
matter power spectrum is computed using the transfer func- 
tion fits given by Eisenstein & Hu (1999). The power spec- 
trum is normalized to the WMAP3 value a% = 0.77 (Spergel 
et al. 2006). 

4. POWER SPECTRA AND MEAN SIGNALS 

Figure 1 shows the y-distortion power spectrum, for differ- 
ent values of the energy injection redshift z !K . The correlated 
term dominates for I < 10 4 , with a broad, relatively flat con- 
tribution between I = 100 and I = 2000, corresponding to an- 
gular scales from 2 degrees down to 5 arcminutes (the angular 
scales on which large scale structure is evident). The Poisson 
term contributes the secondary peak around / = 3 x 10 4 , at an 
angular scale of around 20" (the characteristic separation of 
galaxies). Reducing the energy input redshift z,-„ from 3 to 2.5 
reduces the power spectrum by roughly a factor of 2, with the 
Poisson term and the correlated contribution being affected 
about equally. The dependence on maximum mass M max is 
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FIG. 1 . — The y-distortion power spectrum, for A/ max = 2 x 10 12 Mq and 
Mmin = 1O"M0. The three solid lines from bottom to top are for z,„ = 2.5, 
2.8, and 3. The signal is a combination of the correlation term (which peaks 
at / = 290) and the Poisson term (peaking at / = 3.6 X 10 4 ). Also plotted is 
the primary anisotropy (dotted line) and the noise level per / value (dashed 
line) for an ACT-like model experiment covering 400 square degrees with 1 
arcminute resolution and a pixel noise of 2 /zK. 

realtively weak: the power spectrum amplitude increases only 
by a factor of around 60% if the maximum mass is increased 
by a factor of 5. For comparison, Fig. 1 also shows the pri- 
mary microwave background anisotropy, and the noise per 
I value for a model ACT-like experiment which maps 400 
square degrees at one arcminute resolution and noise of 2 /iK 
per pixel, using the approximate formula in Jungman et al. 
1996b. Note the signal discussed here is above the noise level 
for a range in I. Other signals will also contribute, particularly 
galaxy cluster SZ distortions (Komatsu and Seljak 2002). 

Root-mean-square temperature fluctuations are obtained by 
converting from y-distortion to temperature power spectrum 
at a given frequency, and then convolving with a gaussian 
beam profile. We consider three different beam widths: an 
ACT-like beam of 1 arcminute, and 2 ALMA resolutions of 
15 and 5 arcseconds. The results are shown in Table 1, for 
the power spectrum with z,„ = 3, M max = 2 x 10 12 M Q , and 
M m in = 1O U M0, corresponding to the solid line in Fig. 1. 
Nominal noise targets for ACT are a few /iK per arcminute 
pixel (Kosowsky et al. 2006); the signal calculated here will 
be at the detector noise threshold. For ALMA, la sensitivities 
in continuum brightness are given in Table 2, for a 50-antenna 
configuration, an integration time of 1 hour and a velocity res- 
olution of 1000 km/s. These are calculated using the sensitiv- 
ity formulae given in ALMA memo 393 (Guilloteau 2003). 
With an integration of this length, the SZ signal from quasars 
is also near the noise limit for a compact configuration with 
15 arcsecond resolution; substantially longer integrations will 
achieve raw sensitivity for detection of the signal. 

5. DISCUSSION 

Quasar feedback is an important issue for cosmological 
evolution because of its potential for providing nongravita- 
tional heating and its global impact on structure formation. 
The resulting Sunyaev-Zeldovich distortion of the microwave 
background due to this feedback process is a robust observa- 
tional consequence of this theoretical scenario. Similar sig- 
nals due to supernovas from the first generation of stars have 
previously been investigated (Oh, Cooray, & Kamionkowski 
2003), with an amplitude the same order of magnitude as in 
Fig. 1 . Both of these signals, arising from the superposition of 
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Temperature 
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145 


60 
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0.09 
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15 


2.32 
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15 


0.11 


265 


15 


1.75 


145 
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2.35 


220 
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0.11 


265 
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1.78 



TABLE 1 

Root-mean-square temperature fluctuations at ACT 

FREQUENCIES AND THREE ANGULAR RESOLUTIONS. 



Frequency 


Resolution 


Baseline 


Sensitivity 


(GHz) 


(arcseconds) 


(km) 


( M K) 


145 


15 


0.0284 


2.41 


145 


5 


0.0853 


21.74 


220 


15 


0.0187 


1.76 


220 


5 


0.0562 


15.84 


265 


15 


0.0156 


1.63 


265 


5 


0.0467 


14.63 



TABLE 2 

ALMA CONTINUUM BRIGHTNESS SENSITIVITIES FOR A ONE-HOUR 
OBSERVATION. 

many small-angle SZ distortions, result in a small net ampli- 
tude temperature fluctuation spread broadly in angular power. 
The signal is quite interesting in its own right, and also gives 
a possible noise source for current arcminute-resolution mi- 
crowave experiments. 

APEX and SPT both aim to produce maps of the microwave 
sky with detector noise of around 10 /iK per arcminute sky 
pixel. The SZ distortions from quasars is well below this noise 
level, and should not be an important contributor to the noise 
budget for these experiments. ACT initially plans to cover 
less sky at greater depth, with a pixel noise of a few /iK. If 
ACT eventually attains a nominal map sensitivity of 2 /iK per 
pixel, the quasar SZ signal will be at the noise limit. Also 
note that infrared emission from high-redshift dusty galax- 
ies, as detected by SCUBA (Borys et al. 2003), will provide 
a confusion noise limit somewhere between 1 and 10 /iK at 
ACT and SPT frequencies and resolutions. While individual 
quasars will not be detected by these experiments, stacking 
microwave images centered on known quasars can give a high 
significance detection. For a map with noise of 10 /iK per 1' 
beam, each quasar SZ signal will have a signal-to-noise ratio 
of around 0.2. Data from the Sloan Digital Sky Survey reveal 
around 50 photometrically detected quasars per square degree 
(Richards et al. 2006); a 200 square-degree microwave survey 
with optical imaging will give around 10 4 quasars. Stacking 
map sections centered on N quasars will increase the signal- 
noise ratio by a factor of N 1 / 2 , giving a 20<r detection of the 
quasar SZ distortions. The actual sky area needed could be 
much smaller, since many low-redshift galaxies once hosted 
quasars, and the cooling time of the hot bubbles surrounding 
them is long. Distinguishing this signal from any intrinsic 
emission may require measurements at multiple frequencies 
or increased angular resolution to resolve the sources. 

Interferometric experiments like ALMA will have superior 
spatial resolution and, due to a large number of telescopes, 
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very high potential sensitivity. ALMA's spatial resolution will 
be sufficient to resolve the intergalactic separation and thus 
see individual galaxies, reducing any confusion noise. In a 
high-resolution configuration, ALMA could distinguish be- 
tween the SZ signal from a hot bubble surrounding a galaxy 
and infrared emission coming directly from dust in the galaxy 
or from a dusty torus surrounding a quasar (Kawakatu et al., 
2007). The ALMA design reference science plan (Dishoeck, 
Wootten, & Guilloteau 2003) proposes over one year of to- 
tal observations devoted to galaxies and cosmology, likely in- 
cluding a deep survey in several bands which would be useful 
for detecting the signal considered here. The angular reso- 
lutions of 5 and 15 arcseconds in Table 2 are a good match 
for the angular size of the hot bubbles, which is around 10 
arcseconds. Substantially larger beam sizes will not result in 
increased sensitivity, since more noise will be added from the 
regions outside of the unresolved sources, while smaller beam 
sizes will give a lower raw sensitivity for a given integration 
time. 

The prospect of clear detection of SZ distortions from 
quasar energy output is exciting. Such observations will probe 



models of heating of the intergalactic medium, constrain the 
thermal history of the universe at low redshifts, and provide 
important inputs into models of galaxy and quasar evolution. 
We hope these results will prompt more detailed investigation 
of both the signal and its detectability with ALMA and up- 
coming microwave experiments. 
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